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O ; Abstract 

j^ , Strongfy interacting binary mixtures of superparamagnetic colloidal particles confined to a two- 

dimensional water-air interface are examined by theory, computer simulation and experiment. The 
P^ ■ mixture exhibits a partial clustering in equilibrium: in the voids of the matrix of unclustered big 

o, 

^ O^, particles, the small particles form subclusters with a sponge-like topology which is accompanied by 

a characteristic small-wave vector peak in the small-small structure factor. This partial clustering 
d^ • is a general phenomenon occurring for strongly coupled negatively non-additive mixtures. 
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Self-organization processes in suspensions of colloidal particles are of key relevance for 
the construction of nanodevices with desired mechanical, rheological and optical properties. 
In particular, colloidal aggregates or clusters can be used as supramolecular building blocks, 
motivated by the quest for photonic band-gap materials yj, molecular-sieves J2|, y] and 
micro-filters with desired porosity |Ji]. Colloidal cluster formation is mainly controlled by 
the effective interparticle forces [5|]: strong attractive forces between the colloidal particles 
will drive them towards irreversible coagulation into dense aggregates while repulsive forces 
in a stable suspension keep the particles apart, avoiding clustering. More recently, it was 



predicted that a combination o: 



short-ranged attractive and long-ranged repulsive forces can 



lead to equilibrium clusters (a III Isl for the case of hard colloids, whereas particular classes 
of ultrasoft colloids can show clustering even in the complete absence of attraction [9|. The 



cluster morphology critically depends on the detai 



10 



s of the interaction, being, e.g., spherical 



ll|. The equilibrium clustering has been 



or worm-like with an internal spiral structure 

confirmed in recent experiments |l2l . Il3l . Il4| of one-component weakly charged suspensions 

with nonadsorbing polymer additives. 

In this Letter, we address equilibrium clustering in two- component suspensions of big and 
small colloidal particles. Our motivation is twofold: first, any realistic sample consists of 
several components and hence considering a binary mixture is an important step towards 
a control of novel composite materials. Second, and more fundamentally, mixtures exhibit 
much richer behavior than their one-component counterparts in general. Therefore, there 
is need to explore whether and how the scenario of equilibrium clustering occurs in two- 
component systems. Specifically, we consider here a binary system of superparamagnetic 
colloidal particles that are confined to a planar water-air interface and exposed to an external 
magnetic field perpendicular to the interface. The magnetic field induces a magnetic dipole 
moment on the particles, resulting into an effective repulsion between all parallel oriented 
dipole moments, which scales with the inverse cube of the particle distance. This enables a 
direct comparison between theories and computer simulations based on a pairwise dipole- 
dipole interaction potential and renders these suspensions into an ideal model system (l5l |. 
allowing thereby to additionally confront the theoretical findings with direct experimental 



las been exploited to study 



data. The advantage of the availability of experimental systems 

two-dimensional melting in the one-coniponent system |l6l . Il7 | and for the study of the 

dynamics [l^ and the glass transition []j| of the binary system. Yet, cluster formation 



for this system has not been studied to-date. The system at hand presents a host of novel 
features: it is a two-dimensional system, it has clearly defined interactions between all its 
constituent particles, which take a simple form and are steered by a single and readily tunable 
external field; it is experimentally realizable and allows direct comparison with theory; and, 
finally, it is apt to direct visual observation. 

We find that colloidal clustering is specific to the components in the repulsive mix- 
ture: subclustering is possible where only the small particles form clusters but the big 
ones are unclustered, despite of the fact that they are strongly interacting with the small 
ones. This partial clustering does not occur in the effective one-component colloid systems 
since the structure of the additives simply follows that of the colloids. We trace the ori- 
gins of the phenomenon to the negative nonadditivity J20| of the mixture, which expresses 
the fact that the cross-interaction between unlike species is less repulsive than the sum 
of the interactions between like species. Quantitatively, negative nonadditivity, A < 0, is 
defined as follows. Let Vij{r), i,j = 1,2 be the three interaction potentials between all 
species and introduce the corresponding Barker-Henderson effective hard core diameters 
aij = Jq dr{l — exp[—Vij{r)/kBT]}, with Boltzmann's constant k-Q and the absolute tem- 
perature T. Then, A = 2o"i2 — (o"ii + 0-22). The partial clustering effect should there: 



be more generic to negatively nonadditive binary mixtures of, e.g., charged colloids 



ore 



2l| or 



star-polymers and linear chains 2^. However, partial clustering is absent in strictly addi- 
tive (A = 0) hard sphere mixtures and in colloid-polymer mixtures, which are positively- 
nonadditive (A > 0). The structural signature of this phenomenon is the appearance of an 
additional small-wave-vector pre-peak in the small-small structure factor correlations, which 
does not show up in the big-big structure factor. Based on a topological analysis using the 
Euler characteristic, we show that the small-particles have a nontrivial, sponge-like struc- 
ture. We base our findings on liquid-integral equation theory, computer simulations and on 
video microscopy data of real samples, finding very good agreement between all. 

In our theoretical model, we consider two different species of colloidal particles moving 
freely in a two-dimensional plane, which corresponds to the experiment (see Fig. Q). Each 
component is characterized by its partial density pi and its susceptibility Xi? "^ = 1)2. The 
external magnetic field B standing perpendicular to the plane, induces in each particle 
a magnetic moment nii = XiB, i = 1,2. The particles are superparamagnetic, i.e., the 
magnetic dipole of each species perfectly aligns with the external field. We model all particles 
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FIG. 1: Schematic view of the setup: binary mixture of superparamagnetic colloidal particles at 
an air-water interface in an external magnetic field B perpendicular to the plane. 



as point-like but refer thereafter to the particles having the larger susceptibility as the 
'big' (species 1) and those with smaller susceptibility as 'small' (species 2). The particles 
interact via purely repulsive dipole-dipole pair potentials [3vij{x) = Tij/x^, i,j = 1, 2, where 
P = {k^T)'^ and x denotes the distance between any two particles scaled over the average 
interparticle separation between big particles, x = r^/pi. The interaction strengths Tij are 
given as Tij = PxiXjB^ /(^ii- The system is fully characterized by three quantities: the 
density ratio P2/P1, the susceptibility ratio X2/X1 < 1 and one of the three interaction 
strengths, which we pick to be Tu in what follows. 

We have performed Monte Carlo computer simulations [2^ in a square box with periodic 

boundary conditions. Finite-size effects were carefully checked for by varying the box size 

rom 500 to 2000 particles. Additionally, we solved the binary Ornstein-Zernike relation 

24 1 in two spatial dimensions, using the Rogers- Young closure (RY) 25[. The RY gives 



reliable results in comparison to the simulation for all parameter combinations investigated, 
see representative results in Fig. El Thus, we rely on this closure to calculate the structure 
of the fluid. In Fig. 121 the clustering of the small particles is manifest by an additional 



small- wave- vector peak in the structure factor 5*22 (fc) whereas Sn{k) is devoid of any such 
signature. Hence, the formation of small particles subclusters in absence of big particle 
clusters is witnessed by structural information in reciprocal space. 

The experimental system is fully described by the theoretical model introduced above 



and is explained in detail elsewhere 
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, |l6[ . We use superparamagnetic particles 26[ with 
susceptibilities xi = 6.2 x 10^^^ Am^/T and X2 = 6.6 x lO^^^Am^/T, which are suspended 
in a free-hanging, flat water droplet attached to a top-sealed glass ring (diameter=8 mm, see 
Fig. Q). Stabilization with Sodium Dodecyl Sulfate prevents the particles from aggregation. 
Due to high mass density (pi = 1.3kg/dm^, p2 = 1.5kg/dm^), both types of particles 
are pinned down to the water-air interface by gravity and form an ideal, two-dimensional 
monolayer of binary dipoles. The relatively small gravitational lengths oili = 8 nm and I2 = 
62 nm for big and small particles, respectively, compared to the particle diameters of 4.7 /xm 
and 2.8 /im ensure an almost perfect realization of a 2D-system. The fact that the centers of 
masses of the different particles do not lie on the same plane can be neglected. Furthermore, 
the flatness of the interface can be controlled in the range of less than 1 pm. Inclination 
control of the whole setup guarantees a nearly horizontal alignment of the fiat surface, 
ruling out the occurrence of any density gradients in the sample. We control the interaction 
strengths T^ between the dipoles by applying an external magnetic field perpendicular to 
the interface. The dipole interaction dominates all other interactions in this colloidal system 
|l5l |. The sample was conserved for months and measurements were taken for a duration of 
up to 24 hours. All necessary data were recorded by video microscopy. Typically, about 1300 
particles were observed in a box of 689 /im x 505 /im, with a total amount of 10^ particles 
in the whole sample. The density ratio P2/P1 was varied between 0.67 and 1. Local and 
statistical properties of the sample are gathered at a rate of about 1 frame/second on all 
relevant time and length scales. 

A comparison between the theoretical and experimental structure factors is shown in Fig. 
01 pertaining to susceptibility ratio X2/X1 = 0.1 and corresponding to the real experimental 
situation. The representative results are shown here for density ratio P2/P1 = 0.89 and 
coupling constant Fn = 4.05. Very good agreement between theory and experiment is 
achieved, supporting the modeling of the system by means of dipolar interactions exclusively; 
evidently, all other residual forces in the colloidal suspension are much weaker and can be 
neglected. A prepeak is seen in the small-particles structure factor 6*22 (fc), whose height is 



comparable with that of the second peak due to the small value of X2- Indeed, as the second 
peak arises from pure small-small interactions, its height is suppressed since the latter scale 
as xi- The prepeak would trivially disappear in both limits X2/X1 -^ (ideal small particles) 
and X2/X1 -^ 1 (one-component system). Again, Sii{k) is deprived of a prepeak. 

Further evidence of the subclustering phenomenon is provided by real-space data. In Fig. 
13] we show experimental and simulation snapshots of the system under consideration. We 
observe a cluster structure of the smaller particles that perfectly fits in the voids shaped by 
the big particles. To better quantify the morphology of the thus forming spatial patterns, 
we employ tools from integral geometrv and in particular the Euler characteristic E, which 



describes the topology of the pattern [27| . We have discretized the snapshot by superimpos- 
ing a grid on the picture, surround each particle by circle and gradually increase the (cover) 
radius Rq of the latter. The Euler characteristic was then calculated on the discrete lattice 
J27[. The resulting spatial pattern consists then of a number C of connected components 
and a number H of holes and E = C — H (27 1. 

In Fig.El^a), we show the ii^-curves of the two components of a binary mixture, compared 
to a one-component and an ideal gas system. The underlying snapshots are taken from 
simulation data. The Euler characteristic of the small component tends to zero much more 
slowly than the ideal gas curve. For the ideal gas, the particles are more or less evenly spread 
over the area, so the vacant space between the particles can be filled much faster than in the 
two-component case where the inter-aggregate space between the small particles is larger. 
A comparison between the big component and the one-component system reveals that the 
Euler characteristic of the former becomes less negative and tends to zero much more slowly 
due to the less and larger voids caused by the aggregation of the small component. Fig.Efb) 
refers to the parameter combination of the experimental sample. There, the E vs. i?o-curves 
from simulation and experiment are compared, showing very good agreement. The prepeak 
position in Fig. Efb) points to a typical cluster size Ld/^/pI = 1.26 which agrees well with 
the position of the most negative value of E2 from Fig. Efb) and is consistent with the large 
holes left there when the covering circle radius Rq = Lci/2. The partial clustering is caused 
by the interplay between the successively stronger repulsive interactions ^22, "^12 and Vn, 
which preclude macrophase separation and favor the nesting of small particles in the voids 
of the big ones. 

We have studied equilibrium cluster formation in a two-component colloidal model sys- 



tern. We find partial clustering: only the small particles show a clustering, as witnessed 
by a first peak in their structure factor while the big particles do not exhibit any peak on 
this scale. The cluster structure does not consist of well-separated isolated clusters with a 
well-defined geometrical shape, such as circular or lamellar, which form for interactions with 
a short-range attraction and a long-range repulsion [7|. Rather, it possesses a percolating 
sponge-like topology, as signaled by a considerably negative £^-value for large covering cir- 
cles. We presented results for Xijxi = 0.5 and 0.1 but the subclustering phenomenon holds 
for a broad range of parameter combinations; more extensive results will be published in the 
future. 

The subclustering behavior is more general and not just limited to paramagnetic colloids. 
Similar effects are encountered in binary star polymer mixtures J28[. Our experimentally 
realizable system offers a clean and reproducible example of this phenomenon which has not 
received due attention in two spatial dimensions up to now. The common feature of systems 
exhibiting partial clustering is that they possess a negative non-additivity, i.e., if their cross- 
interaction is less repulsive than the sum of the two direct interactions. Partial clustering, 

n 

however, is absent in additive hard sphere mixtures |22l]. In systems with positive non- 
additivity like the Asakura-Oosawa-Vrij model J30|, partial clustering was never observed, 
due to the attractions induced by the polymers, which drive macro-phase separation. It 
would be interesting to study the viscosity and the propagation of light in partial clustered 
structures in order to extract their specific material properties. It is further tempting to 
shock-freeze the clustered structure and use it as a bicontinuous device of controlled, random 
nano-porosity. 
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FIG. 2: Partial structure factors Sii{k) for the big particles and 5*22 (A;) for the small particles. 
Computer simulation results are compared to the RY closure for the following parameter combi- 
nation: Til = 4.0, P2/P1 = 0.5, X2/X1 = 0.5. 




FIG. 3: Partial structure factors Sii{k) for the big particles and 5*22 (A;) for the small particles. 
Experimental results (EXP) are compared to theoretical ones (RY) for the parameters: Tn = 4.05, 
P2/P1 = 0.89, X2/X1 = 0.1. 
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FIG. 4: Snapshots of the binary magnetic mixture from experiment and simulation for the param- 
eter combination Tu = 4.05, P2/P1 = 0.89, X2/X1 = 0.1. The big particles are denoted gray and 
the small ones black. Clockwise from the upper left panel: experiment, simulation (both species), 
simulation (small particles only) and simulation (big particles only). 
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FIG. 5: (a) Euler characteristic per particle, E/N , from simulation for a binary magnetic mixture 

and Til = 3.0, P2//O1 = 0.5, X2/X1 = 0-5; (b) Comparison of the theoretical (black lines) and 

experimental (gray lines) of the same quantity for the parameters Fn = 4.05, pil P\ = 0.89, 

X2/X1 = 0.1. Solid lines: big particles; dashed lines: small particles. 
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